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Abstract 

We investigate the two-dimensional (2D) highly spin-polarized electron accumulation layers com- 
monly appearing near the surface of n-type polar semiconductors BiTeX (X = I, Br, and CI) by 
angular-resolved photoemission spectroscopy. Due to the polarity and the strong spin-orbit inter- 
action built in the bulk atomic configurations, the quantized conduction-band subbands show giant 
Rashba-type spin-splitting. The characteristic 2D confinement effect is clearly observed also in the 
valence-bands down to the binding energy of 4 eV. The X-dependent Rashba spin-orbit coupling 
is directly estimated from the observed spin-split subbands, which roughly scales with the inverse 
of the band-gap size in BiTeX. 
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In a nonmagnetic solid compound, spatial-inversion and time-reversal symmetries require 
the degeneracy of an electron's spin, which can be lifted in a noncentrosymmetric system 
by introducing spin-orbit interaction (SOI). The Bychkov-Rashba (BR) effect describes such 
SOI in an inversion-symmetry broken system, by simply considering a two-dimensional (2D) 
electron gas with an electric field applied along its perpendicular direction The Hamilto- 
nian of the BR model is given by Hbr = A o-[E z x k), where A is the coupling constant reflect- 
ing the SOI, a and k are the spin and momentum of electrons, and E z is the applied electric 
field. It gives an the effective magnetic field acting on the electron spin a, whose orientation 
is normal to E z and k. The energy dispersion is characterized by the fc-linear splitting of the 
free electron parabolic band dispersion as denoted by E ± (k) = h 2 k 2 /2m* ±a^\k\, where m* 
is the effective mass of the electron and gives the degree of Rashba-type spin-splitting 
i.e. Rashba parameter. Until now, such asymmetry-induced spin-splitting has been actively 
investigated and designed in semiconductor heterostructures such as the InGaAs/GaAs sys- 
tent g More recently, attempts to induce spin-splitting by applying the electric field in 
centrosymmetric heavy-element oxide semiconductors, such as KTa03, are also reported [3[. 
The spin-splittings in these systems have remained small compared to the thermal energy 
at room temperature, nevertheless, due to the weak SOI and/or small asymmetric E z that 
can be applied by such field effects. 

Previously we reported the findingof a giant Rashba-type spin-splitting realized in the 
bulk of a polar semiconductor, BiTel |4|, la, la, LZJ] - Its peculiar polar crystal structure leads to 
a bulk-induced Rashba-type SOI, similarly as originally discussed in wurzite semiconductor 
materials js|. The Rashba parameter or for BiTel is found to be as large as 3.5 eVA, which 
is comparable to the largest value reported thus far for a surface alloyed metal Bi/Ag(lll) 
system {9]. By utilizing such a large Rashba-type spin-splitting nature built in the bulk semi- 
conductor material, one may obtain a 2D electron system under a strong Rashba field near 
its surface and interface, which has been highly desired from the viewpoint of developing 
various spintronic functions 10, QjJ and multi-functional interfaces 12j . In this Letter, we 
investigate the near-surface regions of the extended families of polar semiconductors BiTeX 
(X=I, Br, and CI) by means of high-resolution angular-resolved photoemission spectroscopy 
(ARPES). Using first-principles calculations combined with the Poisson-Schrodinger (PS) 



method 



13], we precisely characterize the 2D confinement effect on electronic structures 



including both conduction- and valence-bands, emerging as a consequence of charge accu- 
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mulation at the surfaces. 

Single crystals of BiTel were grown by the Bridgman method [4|, whereas those of BiTeBr 
and BiTeCl were grown by the chemical vapor transport method. Laser-ARPES was per- 
formed u sing a vacuum ultraviolet laser {hv = 6.994 eV) and a VG-Scienta R4000 electron 
analyzer 14j with a total resolution energy of 4 meV. /^-dependent ARPES measurement 
{hv = 60 - 84 eV) was performed by using the VG-Scienta SES2002 electron analyzer at BL28 
in Photon Factory, KEK, with a total resolution energy of 20 - 30 meV. Samples were cleaved 
in situ at around room temperature and measured at 10-15 K. Electronic-structure calcula- 
tions for BiTeX were carried out within the context of density functional theory using the 
full-potential augmented plane-wave plus local orbital (APW-LO) method 
the effect of band bending, a 60 layer tight-binding supercell was constructed by down 
the APW-LO Hamiltonian, using maximally localized Wannier functions 



4J. To simulate 
biding 
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13| . The 



band bending potential was then obtained by solving the coupled PS equation [13], assuming 
static dielectric constants 15, 10, and 6.5 for X=I, Br, and CI, respectively. 

The crystal structures of BiTeX are shown in Figs. 1(a) (X = I and Br; space group 
P3ml) and 1(b) (X = CI; P&^mc). They are characterized by the alternative stacking of 
Bi-, Te- and X-layers, which gives rise to the polarity along the stacking axis. In contrast 



to the previous report indicating the disorder of Te and Br sites for BiTeBr [18] , the single 
crystalline BiTeBr we obtained was confirmed to possess the well-defined stacking of Te- 
and Br- triangular lattices by x-ray diffraction measurements. 

The calculated band structures for BiTel and BiTeBr (BiTeCl) as shown in Fig. l(f)-l(h) 
are composed of 6 (12) conduction bands above the Fermi level (Ep) and 12 (24) valence 
bands below Ep. The conduction bands for BiTeX are predominantly composed of Bi 6p 
components. The valence bands, on the other hand, show a systematic X-dependence, as 
seen by the lower half of them moving toward deeper binding energy (Eb) upon increasing 
the electronegativity of X from I, Br, to CI. This indicates that the lower set of the valence 
bands is mainly dominated by X ions, whereas the upper ones primarily reflect the Te 5p 
character. The resistivity for each sample, as shown in Fig. 1(e), indicates metallic behavior 
down to 2 K. The electron densities of X=I, Br, and CI samples are estimated from the 
Hall coefficient to be nbuik = 4 x 10 19 , 4 x 10 18 , and 2 x 10 18 cm -3 , respectively, indicating 
that they all behave as n-type degenerate semiconductors due to slight nonstoichiometry. 
From the band calculation, the Fermi levels of the bulk state corresponding to these carrier 
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concentrations should be located at 142 meV (X=I), 27 meV (X=Br), and 5 meV (X=C1) 
above their respective conduction band minimums (CBMs). 

Figure 2(a) shows the laser-ARPES image for BiTeBr, recorded along the T-K (k y ) di- 
rection. The CBM reaches E-q ~ 0.47 eV, with the remarkable Rashba-type band splitting 
of up to ~ 0.2 eV. The observed CBM is quite deep below Ep, considering that the ex- 
pected Fermi level in the bulk BiTeBr should be located at merely ~ 27 meV above the 
CBM at k z = 7r/c. One can also notice that there are two additional sets of bands at 
shallower energies with CBMs at around E-q = 0.2 and 0.1 eV, thus forming a ladder of 
three Rashba-split subbands. These clearly indicate that downward band bending occurs 
near the surface of BiTeBr, forming an n-type accumulation layer. The conduction electrons 
confined in such 2D accumulation layers are known to form quantized conduction subbands 
I3I . [lgj 2o| . In the present case of BiTeBr, all the subbands clearly distinguished indicate 
the giant Rashba-like splittings, thus offering direct evidence of a 2D spin-polarized electron 
system dominated by the huge Rashba-type effective magnetic field. 

To grasp the topology and shape of these spin-split subbands, the constant energy cross- 
sectional views of the band dispersions are shown in Fig. 2(b). At E-g = 0, i.e. Ep, six 
Fermi surfaces (FSs) corresponding to inner and outer FSs for the n = 1 ~ 3 subbands 
are observed. The large ones, namely the outer and inner n = 1 FSs, are strongly distorted 
compared to the other small circular FSs. Such anisotropy arises from the higher-order terms 
of SOI reflecting the underlying crystal symmetry, beyond the Rashba-type (/c-linear) spin- 
splitting, that becomes appreciable at large k. At deeper E-&, the cross-sections of the bands 
become smaller and more isotropic. The band crossing point (Rashba point) is also clearly 
distinguished as a sharp point appearing at the zone center for some specific E^ (e.g. 160 
meV for the n = 2 subband). To discuss the realized 2D electronic structure, we introduce 
the PS method to calculate the subband structures by taking into account the band bending 
potential and 2D confinement at the accumulation layer 13J. Here, we fix the bulk carrier 
density and calculate the depth (^-dependent energy V(z) and electron density n(z), which 
most accurately reproduce the observed subbands configuration [see Fig. 3(b)]. V(z) and 
n(z) are shown in Figs. 2(c) and (d), which indicate that the electron accumulation layer of 
~ 1 nm thickness is formed at the near-surface region of BiTeBr. 

The systematic results of laser-ARPES and corresponding PS simulation for BiTeX are 
shown in Fig. 3. The experimentally obtained subband dispersions along k y directions are 
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shown in Figs. 3(a)-3(c). Markers showing the peak positions of ARPES intensity estimated 
from the energy distribution curves and momentum distribution curves are overlaid on the 
images. BiTel (BiTeBr and BiTeCl) exhibits two (three) sets of Rashba-type spin-split 
conduction subbands. Subbands of n = 1, 2, and 3 are depicted by red, blue, and green 
markers. The CBMs of the deepest subband (n = 1) are located at E-q ~ 0.27, 0.47, and 
0.41 eV for X=I, Br, and CI, respectively, which are all much deeper when compared to the 
respective bulk E F - E C bm (0.142, 0.027, and 0.005 eV). The 2D electron density at the 
accumulation layer can be estimated experimentally by counting the areas of all the observed 
FSs in Figs. 3(g)-3(i); n 2D = 3.0 x 10 13 , 7.1 x 10 13 , and 4.6 x 10 13 cm" 2 for X=I, Br, and CI, 
respectively. The V(z) and n(z) which most well reproduce the observed subband positions 
for BiTeX are given in Figs. 2(c) and 2(d). The calculated subband dispersions in Figs. 
3(d)-3(f) agree well with the ARPES data, including the group velocities and spin-splittings 
of the respective subbands. By comparing the FSs [Figs. 3(g)-3(l)], we can also demonstrate 
that the PS method quantifies their sizes and shapes (k x -k y anisotropy) as well. Note that 
the n = 1 inner and outer FSs get closer at k x (r~M) direction as compared to k y (T- 
K) direction. This behavior can be qualitatively understood by looking at the bulk band 
dispersions near the CBM. For BiTeBr as shown in Fig. 1(g), for example, the spin-split 
band dispersions along the A-H (k y ) direction remain almost parallel, whereas those along 
A-L (k x ) cross each other at ~ 0.8 eV above the CBM due to the trigonal warping effect, that 
dominates at large k vectors 2l\ . This anisotropy gives rise to the observed deformation of 
the FSs. 

To quantitatively evaluate the spin-splittings, we estimate the Rashba-parameter or for 
BiTeX. Here we calculate «r by using the momentum offset at CBMs (fco) and the Rashba 
energy (£r), as depicted in Fig. 3(d). They are related by £r = h 2 ko 2 /2m* and or = 
2Eft/k in the 2D free-electron BR model. Estimated values for fc and «r are shown in 
Table I. There is no clear n-dependence in the obtained «r (e.g. for BiTeBr, «r = 2.0 ±0.7, 
2.0 ± 0.5, and 2.1 ± 0.7 eVA for n — 1, 2, and 3). This indicates that the Rashba-type spin- 
splitting of BiTeX predominantly originates from the bulk crystal structure, not the potential 
gradient arising from the band-bending. Regarding the X-dependence of Rashba parameters, 
the experimentally obtained values of or for BiTel (3.9 - 4.3 eVA) are significantly larger 
when compared to BiTeBr (2.0 - 2.1 eVA) and BiTeCl (1.7 - 2.2 eVA). Similar values and 
tendencies are also indicated in the series of a R estimated from the calculated subbands, 
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as shown in Table I. The huge spin-splitting of the CB in bulk BiTel has been successfully 
explained by means of k-p perturbation theory {J. According to this model, the SOI- induced 
band splitting near the E-p, i.e. CBM and valence band maximum (VBM), is strongly affected 
by the unusual electronic structures around its band gap. In BiTeX, the VBM (CBM) is 
commonly composed of Te 5p (Bi 6p) orbitals. The band-gap size for BiTeX, on the other 
hand, is known to depend on X anions: The minimum gaps realized in BiTeX are found 



to be 0.36, 0.6, and 0.7 eV for X = I, Br, and CI, respectively, from optical studies [22 ] . 
The narrow energy separation between CBM and VBM, combined with the large atomic 
SOI of Bi, enhances the spin splitting via the second order perturbative k ■ p Hamiltonian 
Considering these aspects, the X-dependent spin-splitting in BiTeX can be primarily 
understood to follow the inverse of the band-gap size, as observed. 

Finally, we investigate the effect of the 2D electron accumulation in a wider energy range, 
including the valence bands. The second-derivative ARPES images for BiTel taken along 
k x direction are shown in Figs. 4(a)-(c) for hv = 64, 72, and 80 eV, which correspond to 
k z = 7r/c, 1.5l7r/c, and (2ir/c). The conduction subbands with Eqbm = 0.3 eV (red) 
and 0.2 eV (blue) for n = 1 and 2, similar to the laser-ARPES result in Fig. 3( 
clearly observed. They are commonly observed for hv = 60 - 84 eV, thus indicating their 
2D characters. As for the valence bands, the peak positions of the ARPES images which 
are fc z -dependent (/^-independent) are carefully separated and plotted by green (red and 
blue) markers. The green markers correspond to the bulk band dispersions, which agree 
well with those previously obtained by the bulk-sensitive soft x-ray /izz-dependent ARPES 
study js} . Because of the surface sensitivity, however, the presently observed ARPES image 
tends to be broadened due to the k z summation. On the other hand, the /^-independent 
2D band dispersions can be clearly observed, whose shapes are strikingly different from 
any k z components of the bulk band dispersions. To account for this, we calculated the 
near-surface electronic structure of BiTel in a wider energy-range. Figure 4(d) shows the 
electronic structure formed in the topmost Te-Bi-I trilayer (TL), whereas Fig. 4(e) indicates 
the projection of all TLs into the bulk. Surprisingly, these calculated data show the ladders 
of spin-split quantized subbands (n = 1 and 2) that are formed near the surface (1TL), 
extending in a wide energy region of the valence bands. By comparing with this calculation, 
the red (blue) markers in the ARPES image can be successfully assigned to the n = 1 
(n = 2) subbands of the valence bands. The observed five sets of VB subband structures 
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(the deepest one could not be clearly distinguished experimentally) agree well with the PS 
calculation. It thus indicates that not only conduction-bands but also valence-bands, that 
are deeply bound to ~ 4 eV below Ep, are well separated from the bulk and quantized near 
the surface as a consequence of 2D electronic confinement. 

The present finding of 2D strong Rashba-field electron system, standing separately from 
bulk in a wide-energy range, thus reveals the potential of polar BiTeX for designing various 
interfaces with strong SOI. It is also worth noting that the Rashba parameter and the 
bulk chemical-potential are tunable in this system by chemical substitution and doping. 
The origin and mechanism of the near-surface band-bending for BiTeX, nevertheless, are 
yet unclear. Very recently, several studies report that accumulation (depletion) tends to 



occur at the Te- (X-) terminated surface 23|, |24j, |25|. It will be interesting to elucidate the 



near-surface charge accumulation in polar semiconductors from the viewpoint of the bulk 
(carrier density, band structure, polarity, etc.) and surface/interface properties, for further 
engineering of 2D strong SOI electron systems. 



In summary, we investigated the electronic structure near the surface of polar semicon- 
ductors BiTeX. The quantized subbands with huge Rashba-type splitting were commonly 
observed in BiTeX, indicating the 2D electron accumulation layer with strong SOI robustly 
realized owing to the peculiar bulk structure. The observed X-dependence of the Rashba 
parameter shows the importance of the band-gap size in designing the SOI field near Ep. 
Together with the PS calculation, the 2D quantization was found to occur in a wide range 
of energy, including the valence bands deeply bound at several eV below Ep. 
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FIG. 1: (Color online) Crystal structures of BiTel, BiTeBr (a) and BiTeCl (b). Brillouin zones 
of BiTel, BiTeBr (c) and BiTeCl (d). (e) Temperature-dependent electrical resistivity of BiTeX. 
(f)-(h) Band dispersions of BiTeX obtained from first-principles calculations. 

FIG. 2: (Color online) (a) Laser-ARPES spectral image of BiTeBr along the k y direction, (b) 
Constant energy cross-sectional views of the subband structure at intervals of 80 meV. (c),(d) Depth 
(z)-dependent potential energy V(z) and electron density n(z) of near-surface BiTeX calculated 
by the PS method. 

FIG. 3: (Color online) (a)-(c) Laser-ARPES spectral images of BiTeX. The peak positions of 
ARPES intensity corresponding to n = 1, 2, and 3 subbands are plotted by red, blue, and green 
markers, (d)-(f) Corresponding subband structures obtained by PS method, (g)-(i) Intensity 
mappings at E-p obtained by laser- ARPES, which correspond to the 2D Fermi surfaces near the 
surface, (j)-(l) Corresponding Fermi surfaces calculated by PS method, z-dependent potential 
energy and electron density in obtained are shown in Figs. 2(c) and (d). 

FIG. 4: (Color online) (a)-(c) Second-derivative ARPES images of BiTel along the k x direction at 
k z = 7r/c, 1.5l7r/c, and (2n/c), respectively. The peak positions of ARPES intensity correspond- 
ing to the fe 2 -dependent bulk bands (^-independent subbands for n = 1 and 2) are plotted by 
green (red and blue) markers. (d),(e) Electronic structures for the topmost Te-Bi-I trilayer and the 
projection of the electronic structures from all trilayers into the bulk, obtained by PS calculation. 
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TABLE I: Rashba energy (meV), momentum offset fco (10 _3 A _1 ), experimental and calculated 
Rashba parameters or and (eVA) of the n-th subband dispersions in BiTeX. Numbers in 
parentheses represent the respective range of error. 





BiTel 


BiTeBr 


BiTeCl 


n 


1 2 


12 3 


12 3 


Er 


108(13) 92(13) 


42(10) 37(7) 35(8) 


25(10) 17(6) 20(6) 


ko 


50(10) 47(7) 


43(10) 38(6) 34(8) 


26(8) 20(6) 18(8) 


OR 


4.3(9) 3.9(8) 


2.0(7) 2.0(5) 2.1(7) 


1.9(10) 1.7(8) 2.2(12) 




3.06 3.02 


1.78 1.67 1.84 


1.28 1.29 1.78 
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Fig. 1 M. Sakano et al. 
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Fig. 2 M. Sakano et al. 
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Fig. 4 M. Sakano et al. 
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